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Abstract 

Quantum decoherence has been studied using nuclear magnetic reso- 
nance(NMR). By choosing one qubit to simulate environment, we examine 
the decoherence behavior of two quantum systems: a one qubit system and 
a two qubit system. The experimental results show agreements with the the- 
oretical predictions. Our experiment schemes can be generalized to the case 
that the environment is composed of multiple qubits. 
PACS number(s):03.67 

1. INTRODUCTION 

Quantum decoherence is a purely quantum-mechanical effect through 
which a system loses its coherence behavior by getting entangled with its 
environment degrees of freedom [1][2]. A real quantum system always in- 
teracts with its surrounding environment. The time evolution induced by 
the interaction introduces entanglement between the system and environ- 
ment when the system initially lies in a superposition of states. Quantum 
decoherence can be viewed as the consequence of such entanglement. When 
the state of the system is described by the reduced density matrix through 
tracing over the environment degrees of freedom, decoherence makes the off- 
diagonal matrix elements approach 0, and leaves diagonal ones unaltered. 
Quantum decoherence is thought as a main obstacle for experimental imple- 
mentations of quantum computation and has been widely studied in theory 
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[3]- [5]. Some theoretical and experimental schemes have also been proposed 
to solve the problem of decoherence in order to fully use coherence in quan- 
tum information [6]- [8]. Myat.C.J. et al observed and studied decoherence by 
coupling an ion in a Paul trap to a reservoir that can be controlled. Their 
experimental results tested their theoretical predictions [9] [10]. 

As a useful tool to research the microscopic world through the macro- 
scopic signal, nuclear magnetic resonance can be used to study the funda- 
mental problems in quantum mechanics. Besides some quantum algorithms 
have been implemented on NMR quantum computers, a quantum harmonic 
oscillator was simulated by S. Somaroo et al [11], and the complementar- 
ity was tested by X.Zhu et al [12]. In this paper, we will study quantum 
decoherence using NMR. Although we choose one nucleus to simulate the 
environment, our schemes can be generalized to the case that the environ- 
ment is composed of multiple nuclei. 

2. The decoherence behavior of a one qubit system 

For the simplest case, we choose a sample containing two coupled spin 
1/2 nuclei to implement experiments. The two nuclei are defined as qubit 
1 and qubit 2 which represent the quantum system and environment, re- 
spectively. The experimental scheme is shown in Fig.l. R(0) denotes a 
rotation operation for qubit 2, and H denotes the Walsh-hadamard trans- 
form for qubit 1 [13]. Experiments start with the pseudo-pure initial state 
I !>i I !>2 by using spatial aver aging [14] [15], where | |>« (i=l or 2) de- 
notes the spin 'down' state of qubit i. H and R(9) transform | |>i | |> 2 
into (| t>i -| |>i)(a(0)| T> 2 +6(0)1 l> 2 )/v / 2, where \a\ 2 + \b\ 2 = 1. After 
the control-not operation is applied, the system and environment lie in the 
entangled state described as 
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where a constant factor is ignored. The state of the system is represented by 
the reduced density matrix[16] 
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where the off-diagonal elements describe the coherence behavior. When 9 is 
chosen as different values, entangled states with different amount of entan- 
glement are gotten [17]. It is clear that the coherence of the system depends 
on 9. 

Our experiments use carbon-13 labelled chloroform dissolved in d6-acetone 
and carbon-13 trichloroethylene(TCE) dissolved in d-chloroform as two sam- 
ples. TCE's structure is shown in Fig. 2. 1 H is decoupled during the ex- 
periments. The two 13 C nuclei are assigned as qubit 1 and qubit 2, respec- 
tively. Data are taken at controlled temperature (22°C) with the Bruker 
DRX 500 MHz spectrometer of Beijing Normal University. For chloroform, 
the resonance frequencies v x = 125.76 MHz for 13 C, and u 2 = 500.13 MHz 
for 1 H . The coupling constant J\i is measured to be 215 Hz. For TCE, 
the resonance frequencies v\ = 125.7869009MHz (117ppm in the spectrum), 
^2 = v\ + 903.6Hz. The coupling constant Jyi = 103.1 Hz. 

The pulse sequence [9} 2 X - [f ]£ 2 - - [n] 1 / - ^ - [f ] 2 transforms the 
initial state | |>i | |> 2 into an entangled state. [9] 2 X denotes a rectangular 
pulse along x-axis for qubit 2, and denotes the evolution caused by 
the magnetic field for when pulses are closed. In NMR experiments, the 
entangled state can be described by the following deviation density matrix[18] 
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sm(|) in Equation(l). By trac- 
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ing p pp over the environment degrees of freedom, the state of the system is 
described as 





(4) 



where the dependence of the coherence on 9 can be represented by a sine 
function. In the NMR spectrum, the sum of the amplitudes of the two peaks 
of 13 C is proportional to the off-diagonal element in equation (4). 

Fig. 3a and Fig.3b show the experimental results when the sample is cho- 
sen as chloroform. When 13 C and X H lie in an entangled state, the amplitudes 
of the two NMR peaks of 13 C are almost equal as shown in Fig.3a, where 
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9 = 50.3°. When 6 is chosen as different values, the dependence of the co- 
herence on 9 is shown in Fig. 3b. The experimental data points can be fitted 
as a sine-shaped curve with the period 1.967T. The error is about 2.0%. It 
mainly results from the imperfection of pulses. When 9 = or ir, 13 C and 
X H lie in a maximally entangled state [18]. If we ask only about the state of 
13 C, it lies in a completely mixed-state, in which the coherence completely 
disappears [16]. When 9 = tt jl or 3n/2 , 13 C and 1 H lie in a product state. 
The coherence of 13 C is maximal. We will discuss this case in detail. If the 
pulse sequence [9]1 — [f]^' 2 is applied to the initial state | |>i | |>2, 13 C and 
1 H are still in a product state. The spectrum of 13 C is shown in Fig. 4a, where 
9 = 50.3°. The amplitudes of two peaks varies versus 9 shown in Fig. 4b, in 
which they are marked by 'A' and 'O 'respectively. Nevertheless, their sum, 
which is corresponding to the coherence of the system, is independent of 9 
as shown by the line marled by '*' in Fig. 4b. 

Experimental results of TCE are shown by Fig.5a and Fig.5b. A carbon 
spectrum is shown in Fig. 5a when 9 = 27°. The dependence of the coherence 
on 9 is shown in Fig.5b. The period of the sine-shape curve is 1.837T. The 
error is about 8.5%. Because the difference of the two 13 C nuclei in TCE is 
much less than the difference of 13 C and 1 H in chloroform, the selective pulses 
used in the latter experiment are less perfect than the former experiment. 
The error is larger when the sample is changed to TCE. 

Our scheme can be generalized to the case that the environment is com- 
posed of N(N > 1) qubits. The quantum network is shown in Fig. 6. For 
convenience, we set the rotation operation as the same form R(9). After the 
system gets entangled with qubits in environment one by one, its state is 
represented as 



When N oo, the coherence approach 0, if sin9 ^ ±1. In our generalized 
scheme, the qubits in environments are independent. We call this mode 
"ideal gas" mode. The result above is in agreement with the one which is 
gotten by harmonic oscillator mode in theory [3]. 

3. The decoherence behavior of a two qubit system in an 

entangled state 

If a two qubit system lies in an entangled state, the off-diagonal elements 
in its density matrix also represent the coherence. We choose carbon 13- 
labelled trichloroethylene dissolved in d-chloroform as a sample. The two 
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13 C nuclei, which are defined as qubit 1 and qubit 2, are assigned as the 
system, and the 1 H, which is defined as qubit 3, is assigned as environment. 
The influence of the other nuclei can be ignored. The pseudo-pure initial 
state | |>i | !> 2 is also prepared by spatial averaging [11]. 1 H is decoupled 
until the entangled state of the system is prepared. We choose composite 
pulse decoupling, and NOE can be ignored due to the short decoupling time. 
The following pulse sequence [§]^' 2 — jj^ — [^}l' 2 — jj^ — [f]^ transforms 
I !>i I !>2 into a basis entangled state [18] 

p s (0) = Ill 2 x - - Iff, (6) 

where 1%, Iy and I z (k = 1,2) are respectively the matrices for the three 
components of the angular momentum of the spins. As soon as the system 
lies in the entangled state, the decoupling pulses are closed. The system and 
environment will evolute in the magnetic field. In the rotating frame, the 
Hamitonian is represented as 

H = Au 12 ll + Au 32 I 3 z + 2nJ 12 I l z I 2 z + 2tt J 23 / 2 / 3 + 2tt J 13 /]/ z 3 (7) 

where Ui/2ir(i = 1,2,3) are the resonance frequencies of spins 1,2, and 3, 
Aco»i2 = 0J\ — oo 2 , ACU32 = ou 3 — u 2 , and H is set to 1. We choose u 2 /2n 
as the frequency of the rotating frame. By applying a hard (nonselective) 
pulse [tt]^' 2 ' 3 in the middle of an evolution period t, the chemical shift evolu- 
tion during this period can be refocused [20] [21], so that the system in the 
entangled state evolutes under 

H eff = 27rJ 12 / z 1 / 2 + 2vrJ 23 4 2 / 3 + 27rJ 13 / z 1 / 3 , (8) 

where the last two terms describe the interaction between the system and 
environment. Only these two terms cause the system to evolute, because 
entangled basis states are the eigenstates of [19]. After evolution time t, 
by tracing over the environment degrees of freedom, the state of the system 
can be represented as 
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where tp 13 = 7rJ 13 t, and ip 23 = tt^23^ In NMR experiments, this equation is 
equivalent to the following one 
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The off-diagonal elements are called double quantum coherence, which can 
be observed through a readout pulse [|]^ which transforms the system into 
the state represented as 
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The coherence behavior can be observed in the NMR spectrum. 

Experimental data are also taken at controlled temperature (22°C) with a 
Bruker DRX 500 MHz spectrometer. The coupling constants J\ 2 = 103.1Hz, 
J 23 = 201.3Hz, and J\ 3 = 9.23Hz. X H is again decoupled during recording 
the FID signal in order to simplify the spectrum. From the view of quantum 
mechanics, decoupling also can be thought as the process of tracing over 
the environment degrees of freedom to get the reduced density matrix of 
the system. We will discuss this point in detail in another paper. Fig. 7a 
is the carbon spectrum through the readout pulse [|]^ when the evolution 
time is 3.50ms. Fig. 7b shows that the amplitude of the left peak of CI (see 
Fig. 2) in the spectrum varies versus the evolution time. The experimental 
data points can be fitted as function 5.8cos(27rt/T), where T = 8.72ms. 
In theory, T = 9.50ms. The error is about 8.2%. It mainly results from 
the imperfection of pulses, and the effect of decoherence which cannot be 
controlled. 

If environment is composed of multiple qubits, the interaction between 
the system and environment is represented as 
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Under the evolution induced by ifj, the state of the system is described as 

7V+2 iV+2 

p 12 (t) = [J cos(<p lk + V2k ) - - I\l 2 y J] cos(<p lk + cp 2k ), (13) 

k=3 fc=3 

where <fi k = nJi k t, and ip 2k = Kj 2k t. If iV — > oo, the coherence usually 
approachs to very fast, if t ^ 0. The decoherence of the system is also 
caused by the entanglement between the system and environment. 

4. Conclusion 

In our experiments, we examine quantum decoherence of two systems 
by two experimental schemes. For the one qubit system, the entanglement 
between the system and environment is realized a quantum network. For 
the two qubit system, such entanglement results from the evolution induced 
by the interaction between the system and environment. Although the two 
schemes seem different, the interaction is their common prerequisite. We 
can get the results when the environment is composed of multiple qubits 
by generalizing our schemes. It is unnecessary to consider the state of the 
environment. In fact, we usually cannot describe the state of environment 
clearly. Through the entanglement between the system and environment, the 
environment degrees of freedom are introduced into the state of the system. 
Decoherence is the result of tracing over such degrees of freedom. If real 
environment consists of multiple particles each of which has more than two 
degrees of freedom, our results are still correct. 
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Figure Captions 



1. Quantum network used to study the decoherence behavior of a one 
qubit system. Qubit 1 and qubit 2 are defined as the system and 
environment, respectively. H denotes the Walsh-hadamard transform, 
and R(9) is a rotation for qubit 2. 

2. The structure of trichloroethylene. The two 13 C nuclei are assigned as 
qubit 1 and qubit 2, respectively. 

3. The experimental results when 13 C and 1 H of chloroform lie in entangled 
states. Fig.3a is the carbon spectrum when 9 = 50.3°. Fig. 3b shows 
the dependence of the coherence on 9. The experimental data points 
can be fitted as a sine-shaped curve with the period 1.967T. 

4. The experimental results when 13 C and X H of chloroform lie in product 
states. Fig. 4a is the carbon spectrum when 9 = 50.3°. The amplitudes 
of two peaks varies versus 9 are also shown in Fig.4b, in which they are 
marked by 'A' and 'O', respectively. Their sum, which is correspond- 
ing to the coherence of the system, is independent of 9 as shown by the 
line marked by '*'. 

5. The experimental results when two 13 C nuclei of trichloroethylene lie in 
entangled states. Fig. 5a is the carbon spectrum when 9 = 27°. Fig. 5b 
shows the dependence of the coherence on 9. The experimental data 
points can be fitted as a sine-shaped curve with the period 1.837T. 1 H 
is decoupled during the whole experiments. 

6. The generalized experimental scheme when environment is composed 
of qubit 2, qubit 3, and qubit N+l. Qubit 1 denotes the system. 

7. The experimental results when the two 13 C nuclei of trichloroethylene in 
the entangled state evolute under if e jj(see the text). Fig. 7a is the car- 
bon spectrum through the readout pulse [|]^ when the evolution time is 
3.50ms. Fig. 7b shows that the amplitude of the left peak of Cl(shown 
in Fig. 2) varies versus the evolution time which is denoted as t. The ex- 
perimental data points can be fitted as function 5.8cos(2nt/T), where 
T = 8.72ms. 1 H is decoupled during recording the FID signal. 

[Figure 1 about here.] 
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